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ABSTRACT 


The residual thermal stresses in 4-layer cross-ply tubes are 
studied. The tubes considered had a small radius to wal 1 -thickness 
ratios and so elasticity solutions were used. The residual thermal 
stress problem was considered to be axi symmetric and three elasticity 
solutions were derived and the results compared with the results using 

classical lamination theory. The comparison illustrates the limitations 
of classical lamination theory. The three elasticity solutions derived 
were: plane stress, plane strain, and generalized plane strain, the 

latter being the most realistic. The study shows that residual stresses 
in both the hoop and axial direction can be significant. Stacking 
arrangement effects the residual stress to some extent, as do the 
material properties of the individual lamina. The benefits of hybrid 
construction are briefly discussed. 
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INTRODUCTION 


The use of composite materials in engineering structures has gained 
wide acceptance in the last decade. The success of the NASA space- 
shuttle program has provided the impetus for the development of large 
structures for space applications. In order to increase payload effi- 
ciency, these structures will have to be fabricated from lightweight, 
highly efficient components. For the following reasons these qualities 
are uniquely offered by composite tubes: (1) Tubes are components of 

high structural efficiency, as measured by torsional and bending stiff- 
ness per unit weight. (2) Composite materials are characterized by high 
stiffness and strength per weight ratios. (3) Composite tubes can be 
fabricated with relative ease compared to other structural members (such 
as I or T section beams) using filament winding, braiding, or other 
manufacturing techniques. (4) Composite tubes eliminate the free-edge 
problem, a problem which exists in other structural components- 

Since space structures operate under zero gravity conditions, they 
are subjected to relatively mild mechanical loads However environmen- 
tal considerations, such as extremely low (-250°F) temperatures, are of 
great concern. The space environment will subject the material to a 
temperature difference as great as -500°F degrees (from the curing 
temperature). This is a substantial temperature differential and will 
undoubtedly result in high residual stresses. Therefore, in order to 
achieve optimal structure design, it is of great importance to be able 
to calculate the thermal stress magnitudes in composite tubes. 
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This paper reports the first step in the investigation of residual 
stresses in composite tubes. Specific interest is in the residual 
stresses in a tube operating at -250°F. With the tube having a cure 
temperature of, say, 250°F, this represents a severe temperature 
effect. In reality the effects of mechanical loads must be superposed 
on these stresses. These loads, however, will be light and the primary 
effect in the tube may well be the thermal effect. Thus the problem 
studied here will be the determination of the thermally induced stresses 
in cross-ply tubes. The problem can be considered axisymmetric. 
Single-material tubes as well as tubes with layers of various materials 
are studied. In the following work, cross-ply tubes are studied using 
an elasticity approach. These solutions are compared with results from 
a laminated shell theory based on Donnell's shell approximation. 

This report begins by presenting a review of the literature rele- 
vant to the present study. Elasticity theory and shell theory solutions 
are reviewed. Next, the governing equations used in the elasticity 
solutions are introduced and the appropriate simplifications and assump- 
tions use to derive the solution for one-layer transversely isotropic 
tube are outlined. The procedure for the solution of a multiple layered 
cross-ply tube is then presented. Following this the shell theory 
solution to the problem is discussed. Finally, numerical results and 
conclusions are presented. 



LITERATURE REVIEW 


The residual stresses in orthotropic hollow cylinders can be 
studied using an exact elasticity solution, an approximate approach 
(such as shell theory), or by numerical methods. Since the elasticity 
approach provides the most accurate solution to the problem, it is a 
very useful tool in determining the accuracy of the approximate and 
numerical methods. Unfortunately, elasticity solutions are more complex 
and thus difficult to obtain. 

In the following discussion, literature related to the thermo-mec- 
hanical response of anisotropic hollow cylinders will be reviewed. 
Specific attention will be directed to the analysis of residual stresses 
in composite tubes. In addition, studies relating to mechanical loading 
of tubes will be mentioned. The method used in investigating mechanical 
loads are similar to those employed in thermal problems. Most of the 
work found in the literature relates to hollow composite tube subjected 
to combined mechanical loads. 

The review begins by discussing the elasticity method. Following 
this is a review of the shell theory approximate method. Finally, the 
two methods will be compared. 

Elasticity Methods 

Sherrer [1] used an elasticity solution to determine the state of 
stress in a multi-layer f i 1 ament -wound cylinder subjected to a combined 
axial, torsion, and pressure loads. In his model the fiber and matrix 
are treated as separate entities. Consequently, his final expressions 
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were more complex than if the fiber and matrix had been treated as 
homogeneous orthotropic material, as is commonly done in the analysis of 
composite materials. Pagano et al . [2], used a plane-strain elasticity 
solution to investigate the uniformity of the state of stress in an 
anisotropic helical-wound composite tube subjected to an extensional 
force. In their study the laminae were treated as homogeneous orthotro- 
pic layers. This resulted in considerably simplified expressions. 

Pagano and Whitney [3] applied a modified plane-strain elasticity solu- 
tion, in combination with shell theory, to study the effect of the 
degree of material anisotropy and shell geometry on the state of stress 
across the cylinder thickness under simulated laboratory loadings and 
end constraints. Their approach consisted of combining elasticity 
solution to an infinitely long cylinder and a shell analysis to a tube 
of finite length. At a certain distance from the constrained ends, the 
shell's stress resultants were assumed to be statically equivalent to 
the stresses calculated by the elasticity solution. Since the problem 
is easily solved by shell theory, this type of analysis simplified the 
boundary value problem. Weng [4] used a plane-strain elasticity method 
to solve the thermal problem in hollow anisotropic cylinders. Weng 
examined the thermal stress in grades ATJ and ZTA graphite cylinders 
subjected to thermal gradients in a radial direction. The bulk graphite 
was -treated as transversely isotropic material with the hoop-radial 
plane taken as the plane of material isotropy. This condition consti- 
tutes a special problem in fiber reinforced composite tubes with axial 
fiber orientation. His solution was similar to the one generally 
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employed in the analysis of isotropic cylinders subjected to radial 
temperature gradients (see Timoshenko and Goodier, page 448 [5]). In 
Weng's solution the constitutive relations were generalized to trans- 
versely isotropic materials. 

The above studies demonstrated that elasticity solutions can be 
obtained for certain anisotropic hollow cylinder problems. However, 
these solutions become more complex for angle-ply cylinders. In addi- 
tion, elasticity solutions to thermal problems involve the solution to 
nonhomogeneous ordinary (or partial) differential equations and thus 
adds further complexities to the problem. In some cases, depending on 
cylinder geometry and loading conditions, one can utilize an approximate 
method such as first or higher-order shell theory in order to simplify 
the solution to anisotropic hollow cylinder problems. 

Shell Approximation Methods 

A number of investigators have utilized shell analysis to quantify 
the state of stresses in hollow anisotropic cylinders under combined 
loads. One advantage of a shell theory approach, as opposed to the 
elasticity method, is the relative ease with which boundary value prob- 
lems can be solved. Whitney and Halpin [6] used Donnell's shell approx- 
imation analysis to characterize the response of laminated composites 
tubes subjected to combined thermal and mechanical loads. Pao et al . 

[7, 8] demonstrated that by applying Flugge's higher order shell theory 
to orthotropic laminated tubes, the thermal loading problem could be 
solved. Pao also compared the results obtained with a higher order 
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solution. He found that the Donnell method was fairly accurate if the 
shell radius-to-thickness ratio, R/t, was large. However, error of up 
to 10% occurred in the hoop force and moment resultants when this ratio 
dropped below 10. Stavsky and Smolash [9] used a shell approximation to 
derive the thermoelastic equations for an orthotropic semi -infinite 
hollow cylinder shell. In addition, they also developed the heat 
transfer equations which govern the temperature gradient through a 
shel 1 . 

It appears that a shell approximate method can be useful in the 
solution of circular hollow composite shell problems. However, one 
should be aware of the various factors which can contribute to error in 
the analysis. In the following section we shall discuss some of the 
factors and compare the elasticity and shell approximation solution 
methods. 

Comparison of the Elasticity and Shell Approximate Methods 

Shell geometry is one factor which determines the accuracy of the 
method employed. Whitney et al. [10] showed that shell approximation 
yielded accurate results, as compared to an elasticity solution, for 
R/t > 10 in a laminated tube subjected to combined mechanical loads. It 
has also been found in previous studies [2, 3, 10] that for large values 
of R/t, the state of stresses predicted in each layer of a laminate tube 
is similar to those predicted for a laminated flat coupon of infinite 
width subjected to the equivalent membrane forces (i.e., using classical 
lamination theory). Rizzo and Vicario [11] used a three-dimensional 



7 


finite element analysis to determine the effect of R/t and length-to- 
radius ratios, L/R, on the state of stress in composite tube speci- 
mens. They found that for a boron-epoxy tube the stress was nearly 
uniform away from the grip (for R/t > 25) and hence could be accurately 
characterized by a classical lamination theory formulation. 

The effect of lamina fiber angle, relative to the tube axis, and 
stacking sequence on the state of stress in a thermally loaded composite 
shell was investigated by Stavsky and Smolash [9]. A similar problem 
was addressed by Pagano and Whitney [3] in composite tubes subjected to 
combined mechanical loads. Both studies found that fiber angle and 
lamination sequence strongly effected the state of stress in composite 
shells under combined mechanical and thermal loads. Whitney [12] and 
Zukas [13] have pointed out that in laminated composite tubes transverse 
thermal expansion as compared to in-plane fiber direction thermal 
expansion is relatively large. Consequently, the assumption of zero 
strain in the radial direction employed by Donnell shell theory will 
result in considerable inaccuracy in the stress calculation. Whitney 
presented a modified shell theory which incorporated the effect of 
transverse normal thermal strain. In comparing his modified method with 
an exact elasticity solution, he found a good level of agreement between 
the two for an angle ply tube of R/t = 10.5. However, an analysis of 
the same tube using Donnell shell theory resulted in a considerable 
amount of error. 

From the literature review, it is clear that for values of R/t less 
than 10, an elasticity solution or a shell theory of higher order than 
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Donnell's shell theory must be used. Whitney has presented a higher 
order shell theory. The elasticity theory Whitney referenced was a 
private communication with Pagano and no written reference was avail- 
able. Thus the elasticity solution in the following section is presen- 
ted for the purpose of documenting the elasticity approach to the prob- 
lem. 



GOVERNING EQUATIONS 

In the following sections the solutions used to compute residual 
stresses in composite tubes are outlined. The discussion begins with a 
description of the laminated composite tube geometry. This is followed 
by a presentation of the equations which govern the behaviour of a 
composite tube. Finally, the analytical solutions to the thermomechan- 
ical problem of composite tubes using various elasticity solutions and a 
classical shell theory approach are presented and compared. 

Geometry 

In the following analyses, three sets of orthogonal axis systems 
will be referred to; polar cylindrical (r, e, and x), laminate global 
(x', y', and z'), and lamina principal (1, 2, and 3) coordinates. 

Figures 1 and 2 show the tube geometry and the coordinate sys- 
tems. The tube's inner radius is denoted as r^ , the outer radius is r^, 
and the tube's mean radius, R. Note that 


2 

and R locates the distance of the tubes wall's geometric midsurface 
relative to the tube centerline. 

The cylindrical r-e-x system measures the location of a material 
point within the tube relative to a reference point on the tube's cen- 
terline. The directions of the r-e-x coordinates are col inear with the 

three base vectors e , e„, and e shown in fig. 1. This local ortho- 

r 0 X 
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F1g. 1. Tube geometry and coordinate systems used. 
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gonal base vector system is located in space by the cylindrical coordi- 
nates r, 0, and x. The quantity r is the radial distance of the point 
from the tube's centerline, 6 is the circumferential position of the 
point, and x measures the position of the point along the tube. 

The laminate global coordinates x', y‘ and z' are, except for the 
prime notation, those commonly used in the study of composite material 
analysis [14]. They are an off-axis coordinate system. The laminate 
x',y',z‘ coordinates are fixed at the laminate's midsurface, i.e., at a 
distance R from the cylinder's centerline. It should be noted that 

r = R - z'. 

The lamina's principal coordinates 1, 2 and 3 coincide with the 
material principal directions. The 1 and 2 axes represent the fiber and 
matrix directions, respectively. The 3 axis is oriented normal to the 
lamina plane and is col inear with the r and z' axis. The orientation of 
the lamina's fiber (axis 1) relative to the global x' axis is measured 
by the angle ^ as shown by Fig. 2. The above conventions were chosen in 
order to unify the notations employed by elasticity and shell theories 
while maintaining the common convention used in the area of laminated 
composite analysis. 

ff 

Constitutive Relations i 

The thermoelastic constitutive relations for an orthotropic lamina 
can be written relative to the lamina coordinate system as 
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where and are the lamina's linear coefficients of thermal 

expansion in the 1, 2 and 3 directions, respectively. The quanti- 
ties e^, e 2 » e^, Y 23 > Yi 3 > Yl2 total strains. The quanti- 

ties .... t ^2 3re the stresses. The lamina's moduli in the 1, 2 and 
3 directions are designated by (i = 1,2,3), the shear moduli for the 
i-j planes by (i,j = 1,2,3), and the material's Poissons' ratios are 
given by v. . (i,j = 1,2,3). The following three reciprocal relations 
hold: 




i ,j = 1,2,3, 


( 2 ) 
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The temperature change with respect to some reference temperature is 
designated by at, where at > 0 corresponds to a temperature rise. The 
square material matrix is designated by [S] and is referred to as the 
compliance matrix, i.e.. 


CS] 
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1 


13 


21 


23 


31 


'^32 


( 3 ) 
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Due to eq. 2, the above matrix is symmetric. The lamina stress-strain 
constitutive relations can be written in an inverted form as: 
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where [C] = [S]"^ and [C] is termed 
inverse of [S] gives the components 

the lamina stiffness matrix, 
of [C] as follows: 
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where 


^ 11 ^ 22^33 “ ^ 11^23 " ^ 22^13 


^33^12^23^13* 


In terms of the engineering constants the C^j components may be 
expressed as follows: 


^ _ ^■^'23^32 , ''l2'^''32''l3 „ ''l3‘^'^12''23 

^11 ■ E2E3C ’ ^12 " £^£36 ’ ^13 " £^£26 


'22 


'^13^31 

^1^3^ 


'23 


V23+V21^13 

£i£ 2C 


'33 


^"^12^21 

£^£26 


( 6 ) 


'44 


= G 


23 


'55 


= G 


13 


'66 


12 


and 


,, _ ^■''l2''2r''23'^32"'’31^13"^''21^32'^13 
^ ~ F F F 

tit2t3 


As a consequence of the fact that S^j = Sj^, it can be shown that = 


Since the equilibrium and kinematic equations of elasticity are 
generally written in cylindrical coordinates, it is desirable to express 
the constitutive relations in terms of the stresses and strains referred 
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to cylindrical coordinate. This is achieved by writing the stresses in 
the 1-2-3 system in terms of the stresses in the x‘-y'-z' system, and 
writing the strains in the 1-2-3 system in terms of strains in the x‘- 
y'-z' system. From that point, transformation of stresses and strains 
from the x'-y'-z' system to the r,e,x system is made. 

The relation between the stresses in the 1-2-3 system and stresses 



with m = cose and n = sine. 
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The relations between the total strains in the two systems is 



The thermal strains transform like the total strains, i.e. 
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where « , ... a , , are coefficients of thermal expansion in the x'-y'- 

X ^ y 

z' system. The coefficients of thermal expansion in the x'-y'-z' system 
are given explicitly by inverting eq. 11. The result is 


2 2 2 2 
“x‘ ~ ot^i = n Oj + m “z' “ “3 

“y'z‘ " " “x'z'* “xV " 2mn(a^-ot2). 


( 12 ) 


Substituting the transformed stresses and strains into the consti' 
tutive equations, eqs. 1 and 4, gives 
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(13) 


where 


m = [T2]'^ CS][Tj] , 


(14) 


CS] being referred to as the transformed compliance. The transformed 
compliance matrix is of the form 
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The inverse relation is given by 



with 

[C] = CT^]'^[C][T2] (18) 


The matrix CC] is referred to as the transformed stiffness matrix. It 
is of the form 
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lamina stiffness coefficients as; 
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Finally, the stresses and strains can be transformed to 
the r-9-x system. The familiar tensor transformation is used, i.e.. 


a*. . 


a 


ik^j^°kA • 


( 21 ) 


The quantity a^^p is the direction cosine between the +m axis in the x‘- 
y’-z' system and the +n axis in the r-0-x system. The stresses on the 
left of eq. 21 are the stresses on the x'-y'-z' system while the stres- 
ses on the right are the stresses on the r-0-x system. The following 
table shows the specific values of a^p in the transformation. 
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Using eqs. 21 and 22, 



Because Tgg and 5"^^, 5”^^, are the only nonzero terms in 
the 4th and 5th row and column of the CC] and Cs] matrix, respectively, 
the constitutive equations in the r-9-x system can be written as 
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With eqs. 25 and 26, the constitutive behavior of the tube is 
expressed in terms of the r-0-x coordinate system. 

Equilibrium 

The stress equilibrium equations can be derived if the equilibrium 
of forces which act on a small element isolated from a cylindrical body 
are considered. Summing forces in the x, 8 and r directions results in 
the following stress equilibrium equations: 
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where Fp, Fg and F^^ are the body forces in the r, 9 and x directions. 
In the present work these forces will be zero. 


Strain-Displacement Relations 

The strain-displacement relations in polar coordinates can be 
determined in a variety of ways. These relations are 


3w ^ 
3r* 9 


1 r9v , 1 
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_ 9u 9w 

— -K T 15; 

3r 3x 
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X9 


_ 9v 1 9u 

- -57 7 ^ ’ 


where w, v and u are the displacement components in 

the e , and e polar coordinate directions. This is a notation 
r o X 

commonly used in the area of elasticity. It should be noted that the 
engineering shear strain is being used here. In most general elasticity 
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problems, displacements are a function of x, e and r and can be expres- 
sed as: 


u = u(x,0,r) 

V = v(x,0,r) (29) 

w = w(x,0,r) . 


Compa t ibility Equations 

For completeness, the compatibility equations for the strains are 
recorded. In many cases these equations are automatically satisfied. 
For thoroughness, however, they should always be checked. The compati- 
bility equations are 
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Boundary Conditions 

To complete the statement of the problem the boundary conditions on 
all tube surfaces must be specified. The boundary value problem may be 
categorized into two types. In the first type the displacements u, v 
and w are prescribed on the body's surfaces and the stresses, strains, 
and displacements within the body are to be determined. In the second 
type, stresses (or tractions) are prescribed on the body's surfaces and 
the stresses, strains, and displacements within the body are to be 
determined. More often in elasticity problems one is faced with mixed 
boundary value problem, a problem for which stresses are prescribed for 
part of the body's surface, and displacements on the remainder. For the 
problem at hand, the boundary conditions will be generally specified as 
on the cylindrical surfaces, i.e., @ r = >q and r = r^ 
either or w is specified, and 
either or u is specified, and 

either or v is specified. 

on the cylinder's ends, i.e., @ x = Xg^^j 
either or u is specified, and 
either x^g or v is specified, and 
either x or w is specified. 

A n 

However, it is possible to specify integrated conditions at the ends of 
the tube, i .e. , 
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J a dA = specified 
A ^ 

] TT a dA = Specified 
A 

^ ^ = ^end 

etc. 

Here A is the annular cross-sectional area of the layer. 

For the most part, solutions to the general field equations and the 
specified boundary conditions do not exist. Hence, physical assumptions 
are made which simplify the mathematics. The mathematically simplier 
problem is then solved and the resulting answers interpreted in the 
context of the physical assumptions. The approximations for determining 
the residual thermal stresses in composite tubes will be stated later. 

Multiple Layers 

Laminated composite tubes are most often constructed of multiple 
layers. Each layer is governed by the field equations discussed in the 
previous sections. Since these layers are part of the tube, certain 
conditions should be satisfied at each layer interface. For the layers 
within the laminate, these interface conditions replace boundary condi- 
tions on the cylindrical surfaces. The interface conditions between the 
kth and (k-1)^^ layers are that all tractions acting on the surface are 
continuous and all displacements are continuous. These conditions can 
be expressed as; 
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k-1 

u 


k-1 


w 


w 


.k-1 


a 


k 


r 


r 





T 


k-1 
r0 * 



“ ’'interface* 


At the end of each layer 

either o or u is specified 

A 

either or v is specified @ x = Xg^^. 

either or w is specified. 


As with the single layer, it is possible to specify integrated 
conditions at the end of the multiple layer tube, i.e.. 


/ a dA = specified 
A ^ 


/ rt 

A 


x0 


dA = specified 


0 X 


^end* 


etc. 

Here A is the annular cross-sectional area of the multiple layers. 


Assumptions and Simplifications for Computing Residual Thermal Stresses 
In the following sections we shall discuss the assumptions and 
simplifications used in determining residual thermal stresses in multi- 
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layer cross-ply composite tubes. 

1) Axi symmetric simplification 

If the tube's material properties and temperature distribution are 
independent of 0, then the displacements u, v, and w are not dependent 
on 0. Hence 


u(r,0,x) = u(r,x) 

v(r,0,x) = v(r,x) (31) 

w(r,0,x) = w(r,x) 

In addition 

H 0. (32) 

As a result, the strains will be independent of 0. Through Hooke's Law, 
the stresses will then be independent of 0. Thus 


Op = a^(r,x), Og = cjg(r,x), = o^(r,x) 




= ‘^'p(r.x). Eg = E9(r,x), = £^(r,x) 


(33) 


Tp0 = Tp0(r,x), = ^nCr'.x), Y^g= T^g(r,x) 
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The stress-equilibrium, strain-displacement, and compatibility equations 
simplify to: 


simplified stress equilibrium 


3a 


3r 

3t 


3r 

3t 


- + i (»r - <»e> ■ 0 


re . ,2 _ r, 

+ + — T « = U 

r r0 


(34) 


3X 


3a 

+ + 1 ^ = 0 
3r 3x '' 


- simplified strain-displacement equations. 


_ 3w _ w _ 9u 

^r ■ 3r ’ ^0 r * ®x ■ 3x 


_ 3v V 
"^r0 ~ 3r " r 


_ 3u ^ 3w 
^xr ~ 3r 3x 


( 35 ) 


= ^ 
^X0 " 3X 
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simplified compatibility equations 


a^e a^Y 

j + 1 C>L=o (36. a) 

9r3x 


a^e- 1 9e 
9.1 X 

.2 r 


9x‘ 


9r 


i!!rix 
^ 9x 


= 0 


(36. b) 


a 



2 



9r 



9r 


(36. c) 



a2 

r 

9p9x 




= 0 


(36. d) 


a 

9r 


iT-(r^ex)! 


9p 


1 _ 
9r9x 


('• \e) = “ 


(36. e) 


2 

ae a'^ 

_n (re.) = 0 ( 36 . f) 

9x 9r9x 


2) Cross-ply assumption 

Only cross-ply laminated cylinders are considered in the present 
study. Then within each layer there is no coupling between extensional 
and shear stresses. The cylindrical lamina constitutive relations for a 
layer are reduced to the following: 
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- simplified 


stress-strain relations 




( 38 ) 


In a single orthotropic layer or in a multiple-layer cross-ply tube 
subjected to a spatially uniform temperature distribution, there are no 
hoop displacements induced by thermal effects, i.e.. 
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V = 0. (39) 

To summarize, the field equations which govern the residual thermal 
stress state in the individual layers in cross-ply composite tubes in a 
spatially uniform temperature field are: 

- equilibrium 

(40. a) 
(40. b) 
(40. c) 


3a 


3r 

3t 


L + i (a - a„) = 0 
r ' r 9 


r0 


3t 


3r 

3t 


xr 


3x 

3a 




3r 


3X 


- + — T 


xr 


= 0 


- strain-displacement 


e 


_ 3w 
r ^ ’ 



e 

X 


_ 3u 
“ -57 


Y 


r9 


0 


Y 


xr 


_ 9u . 9w 
" TF ^ 


Y 


X9 


( 41 ) 


0 
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- compatibility 


. 2 ^ . 2 ^ 

d G d e d Y 

r X f'X _ n 

^ rr? ^ 


3x 


3r‘ 


3^e„ . 3 g 

9.1 X 

2 ' ~~ m-i --- 

r 


3r3x 
3y 


1 


1121 = 0 


3x 
3 

3r 3r 


3r 

Ir J - r 


3x 


3e 


= 0 


dr 


0=0 (automatically satisfied) 


0=0 (automatically satisfied) 



(42) 


( 43 ) 


or 
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The above system of governing partial differential equations, and 
associated boundary conditions, are still quite difficult to solve for a 
general axisymmetric situation. The most difficult aspect of the equa- 
tions as they relate to tubes is the ability to deal with finite-length 
geometries, i.e., dealing with boundary conditions at the ends of the 
tube. Further assumptions can be made, however, and particular problems 
can be solved. With a proper selection of particular problems, the 
solution to the more general problem can be approximated and perhaps 
bounded. This is the approach to the problem of interest here. Three 
related and solveable problems are examined with the idea that the 
residual stresses calculated from the three simplier problems closely 
approximate and indeed bound the stresses for a more general situation. 
The three solutions are 

i) plane-strain solution, 

ii) plane-stress solution, 

iii) generalized plane-strain solution. 



37 


For certain stress components, the first two solutions represent bounds 
to the problem. The generalized plane strain solution is the most 
accurate representation of away-f rom-end response. The solutions are 
discussed separately below. 

Plane-Strain Solution 

In the plane-strain analysis it is assumed that the tube is 
restrained from axial motion. Consequently, the axial displacement, u, 
is equal to zero and from the strain-displacement relations 


e 


X 



In addition, it will be assumed that none of the other variables in the 
problem vary with x. As a result, the lamina constitutive relations, 
eq. 44, become 


1 0 , 


^ V 

- ot„AT 1 

\ X 


X I 

% 1 

= 


\ 

1 

1 

r r ' 


^xr ^55^xr‘ 


(45) 


When Og and from eq. 45 are substituted into the first equilibrium 
equation, eq. 40a, the result is 


da , _ dc 

+ * ^23 “ 


de. 


+ Co 


dr 


^ 33 . 

dr dr 


+ ^ {(^ 23 " ^22)^0 (C33- C 23 )e^} = p I . 


( 46 ) 
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where: 

I ~ ^ ^^23~ ^^33~ ^23^**r "*" ^^13” * (^^) 


From the strain-displacement relations, eqs. 41, substitution 

for Eg and in eq. 46 leads to a nonhomogenous ordinary differential 

equation for w(r). That equation is 


d^w , dw _ 

^ 33 [ 2 ^ T ^ ~ ^ 

dr^ ^ dr 


22 



(48) 


The solution to eq. 48 can be obtained by making the following variable 
changes: 



This transforms eq. 48 into the following form: 


d^w 




The above ordinary differential equation has the following homogeneous 
solution: 


X-t X„t 
w^(t) = A^e + A^e 


where. 
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To find the particular solution it is assumed 


w 


P 



This yields the following total solution 


^lt ^pt 4. 

w(t) = Aj^e + A 26 + o'" 


C33- C22 


In the r domain 


X, X I 

w(r) = Aj^r + Apr + r . 


r - r 

^33 ^22 


(49) 


In the case when C 33 = C 22 the last term of eq. 49 should be replaced by 


2Z. 


22 


r Inr, and 2 


+ 1 . 


In addition, if and ~ ^ 13 * in the case of 0° lamina 

which assumed to be transversely isotropic, then eq. 49 takes the fol- 
lowing form 


w(r) = Aj^r + A 2 /r . 


(50) 
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The strains and Gq can be determined from the strain-displacement 
relations. The strains are 


e - dw 
r dr 


hr ^22 ^=1 

I 




^9 = F w = 


^ 33 “ ^22 


2 X„-l 
I 


(51) 


+ I A„ r " 


The stresses can be determined by the constitutive relations, namely 


a . 

1 


I - C..a. AT + (C-2 + X,C,^) A, r^^ ^ , (52) 


C33- C22 


A=1 


'^13- 


where 

i = x,0,r 
i = j = 1,2,3 


for stresses and coefficients of 

thermal expansion, 

for stiffness in the same order. 


and the double subscript j indicates summation from 1 to 3. For a 
transversly isotropic 0° lamina, the strains are 



( 53 ) 


and the stresses are 


=6 “ *1 * 
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i and j having the meanings of eq. 52. 

Since is not considered to be a function of x, the third equi- 
librium equation, eq. 40c, results in 

= Br'^ . (55) 

For a single layer the constants Aj and A£ are determined by the 
condition that = 0 at the tube's inner and outer boundary, r^- and r^, 
respectively. With no traction t applied to the inner or outer bound- 

I /> 

ary, B of eq. 56 must be zero. Thus only radial, houp, and axial stres- 
ses are generated. 


Multiple Layer Plane-Strain Solution 

Since the plane-strain solution constrains u to be zero, continuity 
of u at the interface between the (k-l)st and the kth layers is auto- 
matically satisfied. The axisymmetry of the problem and the uniform 
temperature distribution forces the tangential displacement, v, to be 
zero for each layer. Therefore, the continuity of v is automatically 
satisfied. The continuity on w leads to 


w*^(r) 


k-1, ^ 
w (r) 


(a r 


"'interface 



42 


The stress t is zero throughout since there is no applied traction to 

I A 

cause it to be nonzero. Thus the B in eq. 55 is zero for each layer. 
Continuity of the only nonzero stress, a^, at the interface results in 


Or(r) = 0 


''interface 


Enforcing at the N-1 interfaces of an N layer cylinder 

yields N-1 equations to solve for the N Aj^'s and N A 2 's of eq. 49 or 50, 

k k-1 

whichever is applicable. Enforcing a = a at each of the N-1 

r r 

interfaces gives N-1 more equations from which to solve for the Aj^'s and 
A 2 ‘s. Enforcing the traction-free condition on at the inner and 
outer radii gives two more equations. Thus all conditions are satisfied 
in the multilayer tube and the N Aj's and N A 2 's can be determined. 

From eqs. 52 or 54 the stresses can be computed. 

Plane-Stress Solution 

In the plane-stress elasticity solution attention is focused on a 
thin axial slice of thickness 1 . It is assumed that A < < R and the 
stress vector in the axial direction is equal to zero. In addition, it 
is assumed that all stress components are independent of x. Since the 
axi symmetry assumption leads to the stresses being independent of 0, the 
plane stress assumption limits the stresses, and strains, to be func- 
tions only of r. In addition, except for the axial displacement, the 
displacements are functions only of r. The assumptions employed in 
plane-stress are summarized below. 
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stresses 


== "xr = ^x0 " Ve " ° 
o, = a^(r) 

Og = Og(r) 


- displacements 

u = u(x,r) 
w = w(r) 

V = 0. 

For the individual layers in a cross-ply tube, Hooke's relations for the 
above case may be written in the following reduce stiffness form: 


<r> 

1 

I.J 



^22 ^23 


sym Q 


33 


^0 " “0^^ 


- a_AT 
r r 


( 56 ) 


where the Q^^'s are the reduced stiffness coefficients which can be 

* sj 

written as: 


Cl .C, . 

Q. . = c, . - U 

1J 1J f 
^11 


i.j = 2,3 . 


The strain-displacement relations can be substituted into the above 
relations and the stresses written in terms of the displacement w as 
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fol lows: 

% ■ ^^22 (? - “6‘T) + tT23(^ - «, 4T) 

( 57 ) 

°r ' ^23 ■ “ 0 ^^^ '*' "^33 * 

When Og and are substituted into the first equilibrium equation, eq. 
40a, the result is 

d^w 

^33 TT ^33 r ^ " ^22 ^ " r ^ 

dr r 

where 

I “ f(^23” ^22^ *0 ^ ^^33~ ^23^ *^r^ * 

The solution to eq. 58 can be obtained by the procedure outlined for the 
plane-strain case. The result is: 

I 

w(r) = A,r + Apr + r , (59) 

Q33- ^22 

where 

As for the plane-strain case, if ^22 “ ^33 ®*1* 

be replaced by: 
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I 

— — r In r, and x. = ± 1 
2»33 

In addition if a. = a , as will be the case for 0° lamina if we assume 
0 r ’ 

transversely isotropic condition, then eq. 59 becomes 


Ap 

w(r) = Aj^r + — . 


(60) 


The strains e. and e are obtained by substitution of the displacement w 
9 r 

into the strain-displacement relations to get 


e0 =-w = 


_ dw 
r ■ dr 


" 1. 


Q33- Q22 


'^ 33 " ^22 


" .1 V/ 




(61) 


The residual stress components are determined from the constitutive 
relation to be 




^ 33 " ^22 


I (Q2i - X, Qi3) A^r 




(62) 


A=1 


where 

i = 9,r for stresses and coefficients of 

thermal expansion, 

i = j = 2,3 for stiffness in the same order, 
and the double subscript j indicates summation from 2 to 3. 
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(63) 

(64) 

The boundary conditions on the ends of the cylinder are that all 
stresses are zero there, a condition imposed by the plane stress assump- 
tion. The constants and A 2 are determined so as to satisfy the 
condition that = 0 on the inner and outer radii. Only radial, 

a„, and hoop, a., stresses are generated, 
r 0 

Multiple-Layer Plane-Stress Solutions 

The plane stress solution for lailtiple layers is unrealistic in 
that it actually models concentric tubes which can slide axially rela- 
tive to one another. An axial restraint would require nonzero stresses 
in the direction of the cylinder axis, a condition not allowed. Thus 
the cylinders have no choice but to slide relative to one another. Thus 
uk ^ at the interfaces. Since v = 0 through, v^ = v'^"^ at all 
interfaces. Enforcing = w at N-1 interfaces of a N layer cylinder 

yields N-1 equations to solve for the N A^'s and N A 2 ‘s. The stress 

k k-1 

is zero throughout by the plane stress assumption. Hence rx 

each interface. The stress is zero throughout also. Enforcing 


For the transversly isotropic case, the strains are given by 

A^ 

= A^- ^2 




and the stresses are 


,, + Q,3)+A2(Q,2- q,3)^ 
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k k-1 

Op = leads to N-1 more equations for the N Aj^'s and N A 2 's. 

Using = 0 at the inner and outer radii gives the final two equations 
from which to solve for the N Aj's and N A 2 's. Once all the Aj^'s and 
A 2 ‘s are found, eqs. 62 or 64 can be used to find the stresses in a 
particular layer. As with a single layer only hoop and radial stresses 
are generated. 

Generalized Plane-Strain 

In the generalized plane-strain analysis it is assumed that the 
tube is very long. Interest centers on the portion of the tube away 
from the ends. It is assumed that in this region none of the strains, 
and consequently none of the stresses, vary with x. Because of the 
axi symmetric nature of the problem none of the strains, nor displace- 
ments, vary with 0 either. By using the first two compatibility equa- 
tions, eq. 36, it can be shown ^ constant. Thus for the general- 
ized plane-strain analysis 

= e°, a constant, = eg(r), = e^(r) 

' " ’'re “ ^x<3 " ° 
u = u(r,x) w = w(r) . 

From Hooke's Law 


"^re ~ ''^xe 


0 
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Hooke's relations for the nonzero stresses may be written as follows: 



Substituting the stress components and into the first equilibrium 
equation gives 


.2 

C 

‘"33 ^ . 2 
dr 


r dr^ ■ ^22 2 
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w = — 
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( 66 ) 


where 

I = I i^i2~ ^13 ) 

and 


I = {(C 


23 


^22^ 


«0 " 


(C 


33 


C23) 


^^13” ^12^ “vJ 


It should be noted that eq. 66 is identical to eq. 48 except for the 
nomenclature. Equation 66 can be solved by the procedure outlined 
previously. The resulting solution is 


w(r) = 


A^r 


A2r 


^33~ 


'22 


(67) 



where 


7 
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In the case that C33 = C22 the last term in eq. 67 has to be replaced 
by: 

I 

r In r and p = + 1 

or ^9^ 

In addition if and Cj^2 ~ ^13* will be the case for 0° lamina 

if it assumed to be transversely isotropic, the radial displacement w 
reduces to: 

w(r) = AjT + A2/r (68) 


The strain is obtained by substitution of w(r) into the strain- 
displacement relation to get: 


F - dw _ 
r ~ dr 






X,-l 
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C33- C22 





(69) 
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The stresses components for the layer are determined from the constitu- 
tive relation to get 



50 
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2 

" I 

A=1 


(Ci2^- 




(70) 


where 

i = x,9,r for stresses and coefficients of 

thermal expansion, 

i = j = 1,2,3 for stiffness in the same order, 
and the double subscript j indicates summation from 1 to 3. 


For the transversly isotropic case, 

^r " ^1 ■ “2 
r 


(n) 


S ' ''i * 

r 


and the stresses are 


®i " ^li^° ■ ^1^^2i'^ ^3i^ A2(C2^- 03^) \ . 


(72) 


Since is not a function of x, eq. 40c yields 

= Br’^ . (73) 

The method for determining the constants e°, Aj^, and A 2 for a 
single layer is slightly different for the generalized plane strain case 
than it was for the two other solutions. The conditions = 0 at the 
inner and outer radii still are enforced. However, to determine e°, an 
integrated condition must be used. This condition states that despite 
thermal effects, there is no net axial force in the tube. That condi- 


tion is written as 
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] a dA = 0. 

A 

Here A is the annular cross-sectional area of the tube. Since t =0 

xr 

at the inner and outer boundary, B of eq. 73 is zero. 

Multiple Layer Generalized Plane-Strain Solution 
Since is a constant in each layer, i.e., e°, enforcing u*^ = u*^“^ 
at the interfaces leads to the conclusion that is same in all 
layers. Call this strain e°. It is an unknown which must be solved 
for. Since v = 0 throughout, v*^ = at each interface is automatic- 
ally satisfied. Enforcing w*^ = w*^“^ at the N-1 interfaces of a N layer 
tube leads to N-1 equations from which to find the N Aj^'s, N A 2 's, 

and e°. Since is zero on the inner and outer radii, it will be zero 
rx 

k k-1 

throughout and so is automatically satisfied at each 

rx rx '' 

k k ~*1 

interface. Enforcing = o ^ at the N-1 interfaces provides N-1 more 

equations. Requiring to be zero at the inner and outer radii 

provides 2 more conditions for the 2N+1 unknowns A^'s, A 2 's, and e°. 
Using an integrated condition that 

0^ dA = 0 (74) 

provides the final equation from which to solve for the unknown con- 
stants. Here A is the annular cross-sectional area of all N layers. 

This integral condition is physically interpreted to mean that under 
thermal loading only, there is no net axial force acting on the tube 
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cross-section. Since the stresses are independent of x in this analy- 
sis, the integral condition actually states the absence of a net mechan- 
ical load at every cross-section. In each individual layer, however, 
there are axial stresses, in addition to hoop and radial stresses. 

Details of Multiple Layer Solution 

To illustrate the details of obtaining a multiple layer solution, 
the necessary steps for the generalized plane-strain analysis of a two 
layer cross-ply tube are presented in this section. 

Assume that the inner layer has its fibers running axially (0°) and 
the layer is considered transversely isotropic in the r-e plane. The 
outer layer has its fibers running circumferentially (90°). This lamina 
can be considered transversely isotropic in the r-x plane but that will 
not affect the analysis. 

As stated earlier, for each layer the only displacement explicitly 
involved in the analysis is the radial displacement. However, to apply 
the boundary conditions both the radial and axial stresses are 
involved. For the inner layer, referred to here as layer no. 1, the 
displacement is 


w^(r) = aJ r + A^ /r 


(75a) 


and the radial and axial stresses are 
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CiaG + A^(C23 + C33) + A2(C23“ C33) ^ 


■ ^23“e 


(75b) 


= C^e° + aJ(C^2'*’ (^13) ^2^^12~ *"13^ 1 

• 42«6 • 


(75c) 


For the outer layer, here layer no. 2, the radial displacement is given 

by 


w^(r) = A^r^ + A^r"^ + 


C33- C22 


r, X = 


C22 


'33 


(76a) 


The stresses are 


Ci 3 £° + Aj(C 23+ AC33)r^"^ + 


'2'^23 ^^^33^ 


^33 

■ (^13“x ^23«0 ■*■ ‘^33“r^^’’’ ^ ^ 


(76b) 


C33- C22 


= Cj^e° + A^(C^2'^ ^C^3)r^"^ + A2(C^2“ ^C^3)r“^^'^^^ 




^^ll“x ^12“e ^ ^13“r)^'^ ^ 


(76c) 


^33' ^22 


where the bar over the stiffness coefficients, C^j, designate a 90‘ 
lamina and unbarred coefficients a 0° lamina. 


If the interface radius is denoted as r^, with the inner and outer 
radius being denoted by r^- and r^, the boundary conditions which must be 
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enforced on the cylindrical surfaces are: 


o,hr,) =0 

(77a) 


(77b) 

V^(ro) = 0 

(77c) 

w^(rj) = w^(r^) . 

(77d) 


Across a cross-section, the net axial force being zero, eq. 74, 
translates into 

r r 

j ^ <jy(r) rdr + 2ir j ° o^Cr) rdr = 0 . (77e) 

^ '^1 

Substituting the stresses and displacements given by eqs. 75 and 76 into 
77a-d, carrying out the integration of 77e, and arranging into matrix 
form leads to the system of linear algebraic equations given on the next 
two pages. Solving these for Aj^, A^^, ^ 2 ^ and e°, the stresses 
can be computed in each layer by back substitution into the proper 
equation. 


Laminated Shell Theory (LST) 

The degree of sophistication employed in the analysis of shells 
depends to a great extent on the type of assumptions made regarding the 
changes in shell curvature and twist. In the present work the Donnell 
shell approximation on curvature and twist changes were employed. Under 
such approximations all terms of changes in curvature and twist, except 
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those involving second derivatives of W, are neglected [15]. The above 
assumptions lead to the following form of curvature changes. 


^X0 


ax^ 


a0^ 


1 a^w 


axae 


(78a) 

(78b) 

(78c) 


where R designates the radial distance to the laminated shell mid- 
plane. The above first order approximation gives good results for 
homogeneous, elastic quasi-shal low shells (10 < R/t). It remains to be 
seen if the above assumption on curvature changes will yield agreement 
with a more exact solution such as plane elasticity. The body analyzed 
was axi symmetric and under the influence of constant temperature, there- 
fore all curvature changes reduced to zero. 

X“Xn~Xo“0 (79) 

^X ^0 ^X0 ' ' 


The cross-ply laminated shell constitutive relations are given by the 
following expressions 
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where e e ° and y ^ are the strains at the shell's midsurface, 
X y xy 

membrane force and moment resultants, the N's and M's, are defined 
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The a^j coefficients are the elements of the inverted A^j matrix. 
Ajj and B,jj matrices are defined by. 
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where 
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33 


Since we are only interested in the thermoelastic shell problem. 


In addition, in the present investigation it was assumed that z/R < < 

1. Hence, the z/R term in eqs. 82 and 83 were neglected. The laminate 
midsurface strains are obtained by substituting and into 
eq. 80. The residual stresses in the layer are then calculated 
through the lamina constitutive relations. 
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It should be noted that the M's are not zero and can be calculated by 
eq. 81 since the M^'s are known. 



RESULTS 


In the following section the numerical results of the study are 
presented. The discussion will focus on a set of selected cases that 
illustrate particularly significant findings in the current study. 

The analysis was conducted using T300/934 and GY70/934 graphite- 
epoxy composite systems. The T300 graphite fibers represent a high 
strength-low modulus fiber and the GY70 fibers are low strength-high 
modulus fibers. The material properties or both composite systems are 
given in Appendix A. The curing and service temperatures were taken as 
250 and -250°F, respectively. This yielded a maximum temperature dif- 
ference of 500°F. Results are presented for all three elasticity solu- 
tion methods and the results are compared to the laminated shell theory 
(LST). 

Only 4 layer tubes are considered. The next 12 figures illustrate 
the effects of changing the location through the wall thickness of the 
two T300/934 0° (axial) and the two T300/934 90° (circumferential) 
layers, e.g., [0/90/0/90]^ vs. [90/0/90/0]|. . The stress levels are 
plotted as a function of location through the wall thickness. The ratio 
of the inner radius to the wall thickness were chosen to be 12.5. This 
represented the geometry of tubes being investigated in another phase of 
this overall study. Following this, the next twel ve f i gures illustrate 
the effect of changing the number of 0° and 90° layers in a tube e.g., 
( 0 / 903 )^ vs. (0/90/0/90)^. 

After the various 4-layer stacking arrangements are discussed, the 
results for a 4 layer T300/934 tube with all layers circumferential 
( 9 O 4 ) are illustrated. Following that, the comparison between LST 
stress predictions and the plane strain elasticity solution is illustra- 
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ted. This will illustrate the reason for pursuing the elasticity solu- 
tion in this study. Finally, the effect of material properties on 
residual stress is demonstrated by varying the 0° and 90° laminae mater- 
ial properties in a (0/90/0/90)^ laminated tube. Using two types of 
fibers, T300 and GY70, the benefit of laminate hybridizing on stress 
levels is shown. 

The axial residual stress, a , in a (0/90/0/90). tube operating at 
500°F below its curing temperature is given in Fig. 3. The apparent 
lack of agreement between the plane-strain, generalized plane-strain, 
and LST analyses can be explained by the different assumptions employed 
in each analysis. (The axial stress in the plane stress analysis is 
zero). In the plane-strain analysis it is assumed that the tube is 
restrained at its ends from axial motion (i.e., u=0 all along the 
tube). Therefore, a net axial force is required at the tube ends to 
meet this condition. Consequently, the net axial stress across a tube 
cross-section is not zero as may be seen from the sum of area under the 
plane-strain axial stress cure. On the other hand, in the generalized 
plane-strain analysis it is assumed that the tube is unrestrained at its 
ends. In addition, it was required that the net axial stress across a 
tube cross-section is equal to zero. Away from the ends of a finite 
length tube the generalized plane-strain assumption closely resembles 
the LST analysis in the x direction. Therefore, these solution 
approaches are expected to give close agreement as indeed they do. It 
should also be noted that the axial stress in the 90° laminae (hoop 
fiber) is tensile and of considerable magnitude (“ 12,000 psi ). Such a 
stress can lead to circumferential cracking in those plies. 

Figure 4 illustrates the hoop stress in the (0/90/0/900).). tube. As 
can be seen, the layers with their fiber in the axial direction are 






Fig. 4.. Residual Hoop Stress in a (0/90/0/90)^ Laminated 
Composite Tube 
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under considerable hoop tensile stress (« 12,000 psi). This represents 
tension in the matrix direction at a level that will undoubtedly cause 
cracking parallel to the fibers. Conversely, the layers with their 
fibers in the circumferential direction are under compression in the 
fiber direction. This should cause no problems as far as cracking or 
other damage. There is not much difference in the three elasticity 
solutions. There is a difference between the LST and the three elas- 
ticity solutions. For the case shown in Fig. 4, the elasticity solu- 
tions and the LST solution differ by 40% at the outer radius. That the 
three elasticity solutions are similar is interesting in its own 
right. The closeness of the three solutions demonstrates that the hoop 
stress is independent of assumptions regarding the stress or strain 
state in the axial direction. Recall, in Fig. 3, the plane stress, 
plane strain, and generalized plane strain were quite different. All 
elasticity solutions predict a gradient in the hoop stress in the 90° 
laminae but very little gradient in the 0° lamina. As was shown by the 
elasticity analyses, the hoop stress in the 90° lamina varies as a 
function of r^"^ and where X = plane-strain 

and generalized plane-strain solutions and X = /Q 22 /Q 33 for the plane 
stress analysis. As for the 0° lamina, assuming a transverse isotropic 
condition, the hoop stress will vary as function of 1 /r . 

Since X « 3.7 in the T300/934 material, the hoop stress in the 90° 
lamina will vary with r much more rapidly than in the 0 ° lamina. 

Figure 5 shows the radial stress distribution in the (0/90/0/90)^ 
tube. The stress is everywhere tension but the magnitude is quite low 
compared to the hoop stress or the through-the-thickness strength of a 
lamina. Again, the three elasticity solutions are quite close to one 


another. 
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Figures 6-14 show other stacking arrangements with two 0° layers 
and two 90° layers. It appears from these figures that the maximum 
axial and hoop tensile stress is unaltered by stacking arrangement. The 
peak compressive stress in the hoop direction does vary with stacking 
arrangement as does the radial stress distribution. Clustering the 0°'s 
and/or 90°'s, Figs. 9-14, does little to alter the peak tensile stress 
in the matrix direction. 

Changing one of the 90° lamina to be a 0° lamina doesn't alter the 
tensile axial stress in the single 90° lamina, Fig. 15, but it does 
significantly increase the compressive hoop stress in that 90° lamina. 
Fig. 16. The tensile matrix stress in the 0° lamina is still about 
12,000 psi. For a { 90 / 03 )^ tube, the radial stress, as shown in Fig. 

17, is everywhere compressive. Repositioning the single 0° lamina. 

Figs. 18-20, has little effect on the peak axial or hoop stresses. 

The effects of having more 90° lamina than 0 ° lamina are shown in 
Figs. 21-26. Still the peak hoop tensile stress in the matrix direction 
is about 12,000 psi. The peak compressive stress in the fiber direction 
is greater than in other cases. It is interesting to note that LST 
predicts the wrong sign for the hoop stress at the inner region of the 
cluster of 90° lamina. 

The residual stresses in a (90^)^ laminated tube are presented in 
Figs. 27 and 28 as can be seen these stresses are not zero as might 
intuitively be expected. The residual stresses arises due to thermal 
expansion mismatch in the radial and hoop directions. For the ( 904 )^ 
case, the plane stress and generalized plane-strain solution predict no 
axial stress. However, the plane strain solution predicts an axial 



67 


stress which is practically constant with r and has a value of « -205U0 
psi . 

In order to demonstrate the effect of tube aspect ratio (R/t) on 
residual stress magnitude, the maximum hoop and radial stresses are 
plotted versus R/t ratios in Fig. 29 and 30 for a (0,90,0,90)^ tube. In 
Fig. 29 the maximum hoop stresses in the 0 and 90 degree laminae calcu- 
lated by generalized plane-strain elasticity are normalized by the value 
determined from laminated shell theory. As indicated by the plot, only 
at very large aspect ratios does LSI give close agreement with an 
elasticity analysis. This is particularly true for the 90 degree 
lamina. The residual radial stress approaches a small value (10 psi) as 
the tube aspect ratio (R/t) increases. However, it should be noted that 
even at small aspect ratios the residual radial stress is less than one 
percent of the hoop stress. 

Finally, Fig. 31 demonstrates the effect of laminated tube hybridi- 
zation. It appears that the use of such hybridization can be beneficial 
in hoop stress reduction. Specifically, the tensile hoop stress in the 
0° plies, of a 4-ply GY70/934 tube, was reduced approximately 20 percent 
by substituting T300 fibers for the 90° laminae. The above process 
maintains high axial stiffness while reducing the tensile hoop stress. 

An opposite effect is noted when the GY70 fibers in the 0° laminae are 
replaced by T300 fibers (i.e., increasing the hoop stress in the 0° 
laminae). The above phenomenon demonstrates the usefulness of hybrid 
composite tubes in order to maximize longitudinal stiffness and minimize 
residual hoop stress. 
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Fig. 8. Residual Radial Stress in a (90/0/90/0). Laminated 
Composite Tube ^ 
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Fig. 10. Residual Hoop Stress in a ( 02 / 902 )^ Laminated 
Composite Tube 
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Fig. 12. Residual Axial Stress in a (0/90p/0)^ Laminated 
Composite Tube ^ 
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Fig, 13 . Residual Hoop Stress in a (0/902/0)^ Laminated 
Composite Tube 
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Fig. 16. Residual Hoop Stress »in a (90/0.,). Laminated 
Composite Tube ^ ^ 
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AT = -500 °F 



Fig. 17 . Residual Radial Stress in a (90/0^). Laminated 
Composite Tube ^ ^ 
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AT = -500 °F 


R/t = 12.5 



Fig. 19. Residual Hoop Stress in a (0/90/0^) Laminated 
Composite Tube ^ 
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Fig. 20. Residual Radial Stress in a (0/90/0p)^ Laminated 
Composite Tube 
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Fig. 22. Residual Hoop Stress in a (0/902)^ Laminated 
Composite Tube 
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Fig. 25. Residual Hoop Stress in a (90o/0). Laminated 
Composite Tube o t. 
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Fig. 29 . Normalized Residual Hoop Stress in a (0/90/0/90)^ 
Laminated Composite Tube vs. R/t Ratios. 
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Fig. 39, Residual Radial Stress in a (0/90/0/90). Laminated 
Composite Tube vs. R/t Ratios. ^ 




CONCLUSION 


In the present section the conclusions that can be drawn from the 

investigation are summarized below. 

(1) It was shown that the laminated composite tube thermomechanics 
problem can be treated accurately by laminated shell theory provi- 
ded that the tube's R/t ratio is relatively large (R/t > 10). 
However, at smaller R/t ratios considerable deviation is observed 
between LSI and plane elasticity solutions, particularly in the 
hoop stress for a 90° lamina. 

(2) The residual radial stress was shown to be negligible as compared 
with the other stress components (i.e. , hoop and axial). 

(3) The data presented allows one to conclude that varying the stacking 
sequence while keeping the number of 0° and 90° laminae constant 
has little effect on the level of the axial stress in these lay- 
ers. The radial stresses are affected. However, varying the 
number of 0° and 90° laminae in the laminated tube will effect both 
these stresses. 

(4) It was shown that tube hybridization may contribute to both a 
decrease or increase in the hoop stress. This effect depends upon 
the placement of particular fibers in given directions. 

(5) All plane elasticity solutions agreed well with each other as far 
as the calculated hoop and radial stresses. However, the plane- 
strain and generalized plane-strain show the best agreement with 
each other. 

(6) The generalized plane-strain and LST analyses showed close agree- 
ment on the residual axial stress. However, the plane-strain 
analysis failed to give similar results. This was attributed to 

the kinematic constraints assumed in this analysis. 
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APPENDIX A 


The following material properties of T300/934 and GY70/934 graph- 
ite-epoxy composites were used in the study. 


T300/934 graphite-epoxy 
Moduli (Msi) 

El = 21.3 
E 2 = 1.32 
~ 0.67 
Poisson's Ratios 

V12 = 

V 23 = 0.49 
thermal expansion 


= - 0.04E-6 in. /in. °F 
<X 2 = 18.70E-6 in. /in. °F 


GY70/934 graphite-epoxy 
moduli (Msi ) 

Ell = 


-22 


0.975 


Gj2 = 0-71 
Poisson's Ratios 


V 12 = v ^3 = 0.129 
V 23 ~ 0.49 


Thermal Expansion 


= -.578E-6 in./in.°F 
a 2 = 14.42E-6 in./in.°F 
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